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Abstract 

The difference between the timelike and spacelike meson form factors is analysed 
in the framework of perturbative QCD with Sudakov effects included. It is found 
that integrable singularities appear but that the asymptotic behavior is the same in 
the timelike and spacelike regions. The approach to asymptotia is quite slow and a 
rather constant enhancement of the timelike value is expected at measurable large 
Q 2 . This is in agreement with the trend shown by experimental data. 
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1 Introduction 



There is now a long history of continuous progress in the understanding of electromag- 
netic form factors at large momentum transfer. After the pioneering works |I] leading to 
the celebrated quark counting rules, the understanding of hard scattering exclusive pro- 
cesses has been solidly founded by Brodsky and Lepage [0. 

A perturbative QCD subprocess scaling like as(Q 2 )/Q 2 in the simplest case of the 
meson form factor is factorized from a wave function-like distribution amplitude 



rQ 

<p(x,Q ) = / ij(x,k T )dk T 



(x being the light cone fraction of momentum carried by the valence quark), the Q 2 
dependence of which is analysed in the renormalization group approach. Although an 
asymptotic expression emerges from this analysis for the x dependence of the distribution: 

lf as OC X(l - X) 

in the meson case, it was quickly understood that the evolution to the asymptotic Q 2 is 
very slow and that indeed some non pertubative input is required to get reliable estimates 
of this distribution amplitude at measurable Q 2 . Thanks to the QCD sum rule approach, 
such a function was proposed by Chernyak and Zhitnitsky ||, which were followed by 
other model dependent proposals (|, |[. 

These developments helped theoretical estimates to get closer to real experimental 
data but a severe criticism |J remarked that most of the contributions to the form factor 
were coming from end-point regions in the x integration, especially when very asymetric 
distribution amplitudes such as those of were used. This is not welcome since one may 
doubt the validity of the perturbative calculation in these regions. The recent work of Li 
and Sterman |7j] solves this problem by proposing a modified factorization formula which 
takes into account Sudakov suppression of elastic scattering for soft gluon exchange. This 
inclusion leads to an enlargement of the domain of applicability of (improved) peturbative 
QCD calculations of exclusive processes. 

The case of timelike form factors has not been much studied theoretically [§, [J . Exper- 
imental data on the proton magnetic form factor Gm{Q 2 ) [10 show a definite difference 



between the spacelike and timelike values at the highest measured Q 2 . A recent analysis 
of the ip nil decay |ll] leads to a similar problem for the pion form factor at O(10GeV 2 ) 



transfer. Note, however, that the experimental extraction of the spacelike form factor has 



been recently suspected to suffer from large uncertainties [ pT2| . 

In this paper, we carefully analyse in the Li-Sterman framework Q] the ratio between 
high Q 2 timelike and spacelike meson form factors. Not surprinsingly, we find that this 
ratio goes asymptotically to 1 but we show that this approach to asymptotia is slow and 
that factors of the order of 2 follow at measurable Q 2 from reasonable assumptions on 
wave functions. 
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2 The spacelike form factor 



In this section, we review the formalism as it has been developed for the spacelike 
case. 



2.1 Hard scattering picture 

The spacelike form factor measures the ability of a pion to absorb a virtual photon 
(carrying a momentum q with q 2 = —Q 2 < 0) while remaining intact. It is defined by the 
formula: 

< ir(j/)\J*\*(p) >= e^F(Q 2 ).(p + p'y, (1) 

where e v is the pion electric charge and momenta are defined in Fig. [I]. 

In the hard scattering regime, that is when Q is very high with respect to the low 
energy scales of the theory (the QCD scale A and the pion mass), Brodsky and Lepage 
have motivated the following three step picture for the process, valid in the light-front 
formalism [jT^, ^ : 

• the pion exhibits a valence quark- ant iquark "soft" (see below) state, 

• which interacts with the hard photon leading to another soft state, 

• which forms the final pion. 

This leads to the convolution formula: 

F(Q 2 ) = ^ * T H * ^out (2) 

and the graphical representation of Fig. |2|. 

The most important feature of this picture is that it separates hard from soft dynamics. 
The amplitude Tjj, the interaction, reflects the hard transformation due to the absorption 
of the photon and is hopefully calculable in perturbative QCD, because the effective 
couplings are small in this regime due to the asymptotic freedom. The amplitude i/j, 
the wave function, which depends on low energy dynamics is outside of the domain of 
applicability of perturbative QCD and is, at present, far from being fully understood 
from the theory. It is however process independent and contains much information on 



confinement dynamics. Factorization proofs legitimate this picture ITJ] 



2.2 Infrared corrections 

The need of a careful factorization is due to the infrared behavior of QCD: techni- 
cally, large logarithms (~ ln(Q/X)) appear in the renormalized one-loop corrections to 
naive "tree-graph" (A is some infrared cut off needed to regularize soft and/or collinear 
divergences). As in the renormalization procedure, if factorization holds, these large cor- 
rections should be absorbed, here in the re-definition of the wave function. The proof 
of factorization and its consequences upon the wave functions are studied in the pattern 
of the renormalization group. Without entering into a detailed discussion, let us sketch 
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the procedure (see Jl5[ for more on this leading logarithms calculation and also for the 
renormalization group treatment). 

The first step is to compute the naive hard amplitude, that is consider the tree graph 
of Fig. |3| , and the three other graphs related to it by C and T symmetries. One finds, 
with notations explained on Fig. 0: 

T H = W™ S C FxQ2 xyQ2 + { l_ l)2 _ l£ i (3) 

where all quark momentum components are kept. Note that we have done the usual 
projection onto the pion S wave state: 4' 7T (p) oc ^7 5 ^ and used the C symmetry of the 
wave function. Cp = 4/3 is the color factor, while as is the QCD effective coupling at 
the renormalization point /x. 

To examine one loop corrections to Tjj, the relevant graphs to consider in light-cone 
gauge are those of Fig. 0J. 

They directly lead to the wave function correction, in the "double logarithms" or 
Sudakov region (namely: A <C |q| <C <C u an d Q being respectively the light- 
cone fraction and transverse gluon momentum relatively to the pion): 

^ k > - % f "^W) - -> k+ <* - ^ k » 

+ Cj_ mu* <fq du + m _ k + _ 

where x = 1 — x and the first term in the difference comes from vertex-like corrections 
and the second one from self energy ones; in the infrared region some partial cancellations 
occur between these corrections, but the cancellation is not complete. 

To pursue this analysis, it is convenient to define the Fourier transform in the trans- 
verse plane: 

^(x,b) = / d 2 ke lkh ^{x,k), (5) 



and to separate transverse and longitudinal variations of the wave function. One finds, 
omitting for the moment the second term in Eq. (El): 



^ (1) (*,b) = £l(J ^a s(q 2 )(e-- b - 1)/ ^(x,b) 
+ % I ^(qV tqb /^(^ (0) (^-«,b)-^)(x,b)). (6) 

This equation contains the typical corrections one has to consider in a hard process 
when dealing with either a big (^> 1/Q) or a small (~ 1/Q) neutral object. 
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2.3 Transverse behavior at large distance 

The transverse behavior at large distance is driven by the first term of the previous 
equation, thanks to the vanishing of the summation with the oscillating components. This 
occurs when 6 = |b| is greater than at least a few times the inverse of the upper bound of 
the corresponding integral: xQ/^/2. As a consequence, in the remaining expression, the 
infrared cut-off A can be replaced by the natural one 1/6, above which the vextex and self 
energy corrections do not compensate one another. Thus we get: 



i(i) = - s(l nw(») xg/ In 1/6 \ 



(7) 



with (3 = (11 — ^p)/4, rif being the number of quark flavors. Here and in the following, 
it is understood that the energies and inverse separations are in the natural Aqcd unit. 

We have kept the single log term which occurs in the integration, because it is the one 
necessary to express the true dominant large 6 suppression, which one obtains in a more 
complete treatment (that is leading and next to leading one) 

After the ressummation of the ladder structure to all order, the above Sudakov factor 
exponentiates. Taking into account the term obtained with the substitution x — > x = l—x, 
we get: 

*P(x, 6, Q) = e -**fi,Q)-«*m$Wfa 6, Q). (8) 

Thus we get a strong suppression of the effective wave function as b — > 1/A, whatever 
the fraction x is, provided that Q is reasonably large. 

The remaining object i})^ is a soft component to start with. It is soft in the sense that 
it does not include loop-corrections harder than 1/6. One may modelize it by including 
some 6 behavior or simply relate it to the distribution amplitude |y| setting: 



ij {0 \x,b)^ [ 1/b tlj(x,k)dk = <p(x;l/b). (9) 



o 



2.4 Transverse behavior at small distances 

The first term in Eq. (^j) is negligible when the oscillating term remains close to 1 
in the range of integration. This happens for 6 a few times less than max _1 (xQ, xQ). In 
this case, soft divergences cancel one another and one finds: 

^W(s) =£— / dx'l- — — — - — —6(x-x') + - — — — - — —9(x'-x)\, (10) 
2 Jo [ x — x' x' — x J 

with the notation: 

1 [Qd\ 2 1 (lnQ\ 

We displayed this equation in a slightly different form than in the large 6 case to 
explicitly show that Eq. (|j), in the limit of small 6, is related to the distribution evolution 
proposed in [ |13|j . 
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Let us shortly review how this comes. Brodsky and Lepage had proposed a simpler 
factorization formula for exclusive processes 0. It is easily derived from the previous 
treatment if one assumes that neither the wave function nor the hard amplitude give im- 
portant contribution to the form factor when the transverse momenta are big. Neglecting 
all transverse momenta in T# leads therefore to consider the /cy-integrated quantity: 

tp{x) = J d 2 kip(x,k) (12) 

This distribution amplitude related to the wave function at b = has a dependence 
in Q associated with the remaining collinear divergences in Eq. fllPl). Indeed, the expo- 
nentiated form of this convolution equation, once it is written for the distribution ip and 
generalized to other regions than the Sudakov one, leads to the celebrated expansion of 
ip(x,Q)/xx in a linear combination of a running logarithm together with a Gegenbauer 
polynomial. However, whereas this slow evolution is predictible, the expansion at some 
finite Q is inaccessible from perturbative reasoning. 



3 The timelike form factor 

3.1 Quark and gluon poles 

In the timelike region, the hard process ruling 7* — > 7r + 7r~ is drawn in Fig. [5] and 
the hard amplitude is simple to deduce from the spacelike formula changing p — > — p" 
or Q 2 — > —W 2 , W 2 = q 2 . The new feature with respect to the spacelike form factor is 
that the contour of transverse momenta integration now goes near poles located at either: 
k 2 = xW 2 + is or: (k — l) 2 = xyW 2 + is. 

These poles are automatically ignored in the pattern of the Brodsky Lepage formal- 
ism as being to far from the contributing region of integration. However, whereas this 
argument is reasonable at asymptotic regime, we can expect some consequences of the 
presence of these singularities when the energy is not so high. 

Technically, these poles are, except in the end point regions (x,y — > 0), far from 
the bounds of integration of the two independant variables k = |k| and K = |k — 1|. 
Therefore, we may evaluate the integral by deforming the contour of integration in the 
complex plane of each of these variables. 

Another question one may worry about, is the physical origin of these poles. A com- 
plete physical amplitude, for example the form factor F(Q 2 ), considered for complex value 
of Q 2 , has poles and cut along the real negative axis reflecting the existence of interme- 
diate physical (on mass shell) states. These intermediate states are hadronic ones and 
therefore correspond to the "asymptotical" objects of confined QCD. The poles we en- 
counter in our present computation, internal gluon or quark lines going on mass shell, of 
course, do not correspond to observable states. However they only appear in a differential 
amplitude which itself is not observable. Provided, this differential amplitude is integrable, 
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the resulting form factor will only contain, as a remainder of this kind of singularities, a 
real and imaginary parts which one would also expect in a purely hadronic computation. 



3.2 Hard scattering amplitude in b space 

The Sudakov suppression is likely to be important at timelike transfer, so it is nec- 
essary to take the transverse Fourier transform of the hard amplitude. Furthermore, 
whenever this is possible, it is interesting to get an expression not only for negative 
t or positive s, but also for complex values of the generalized transfer. Let us define 
z = \f—i , arg(z) G — | , so that in the spacelike side of the complex plane, we get: 
z = Q whereas in the timelike side: z = —iW. 

The expression of the form factor, with the Fourier transform of Eq. ([3]) and the 
replacement of Q by z, is: 

F = 16itC f J dxdy J b 1 db 1 ^{x,b 1 )b 2 db 2 ip{y,b 2 ) a s T(b 1 ,b 2 , x,y, z), 

T = K (^zb 2 ) xz 2 {9(bx - b 2 )I (^zb 2 )K (^zbx) + (b 2 <-> 6i)} , (13) 

where angular integrations have been done thanks to the cylindrical symmetry of both 
hard amplitude and S wave wave function. b 2 and \b\ — b 2 \ are the transverse distances 
of, respectively, the gluon vertex and the internal quark vertex. The functions Kq and Iq 
are modified Bessel functions of order 0, the first appearing in the equation comes from 
the gluon propagator, while the remainder comes from the quark propagator. 



3.3 Asymptotic behavior 

Bessel functions have different asymptotic behaviors in various directions of the com- 
plex plane: for \(\ — > oo, arg(£) G 



-- 

2 u 



we have HHJ 



*o(C) « xh?e 



^e-^o(C)-gcosh( C + ^), 



(14) 



we thus have to study how the asymptotic dependence of the form factor is affected by 
this direction. In particular in the timelike limit, the integrand is no more exponentially 
suppressed. 

There is, a priori, no general constraint to ensure that the limit of some observable 
like a form factor should be the same in every directions in the complex plane. Even 
though F(z) is analytic and: 



F{z) oc 



a s {zr 



(15) 



with the limit: e(z) — »j^t^=o° 0, our ignorance of the true form of e prevents us from 
concluding when another direction is considered. However, because one expects that 
the same kind of physics underlies exclusive processes, we expect that, at least in an 
asymptotic regime, we should find, for the leading behavior, the overlapping of exactly 
the same soft and hard amplitudes. 
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As a first step in the understanding of the features of the whole form factor, we may 
evaluate the b integrals in an analytical way, putting a simple form for the wave function: 

#(b)=9(B-b), (16) 

which automatically provides a cut-off to avoid indefinite integral. We also forget here 



the possible running of the coupling with transverse distances (see subsection |3.4| ) which 
appears with Sudakov suppression, here ignored. 

With these simplifications, one finds for the integral over 61 and b 2 '■ 



bidbiijj (^^2^2^ (b 2 )T(b u b 2} x, y, z) 

B 



xy z 2 z 3 / 2 



f(x,y,zB), (17) 



with / a function that we refrain from quoting here due to its lack of interest, except for 
its generic behavior for large xy\z\B: f ~ e~ zB . 

As long as we avoid the timelike limit, we find the following leading behavior for the 
integration in the transverse plane: 

1 - - 



xy z 2 



which displays the expected selection of small configuration by the hard process. 

In the timelike region, this is no more true, as we get a modified power dependence: 



Be", (19) 



with the appearance of the size B, together with an oscillating factor. Of course, we 
suspect here that we have found essentially the limit of our model object; nevertheless, 
we may anticipate that some reminiscence of this rather different behavior will occur in 
the non asymptotic regime. 

A source of modification to the above result is the transverse behavior of the wave 
function. The rectangular form which we have choosen above and its steep variation 
reduces the occurence of cancellations expected with an oscillating integrand. If we were 
speaking of Fourier-transform we would say that a rectangular function has relatively large 
components at large momenta in comparison with any similar but smoother function. 

The Sudakov behavior reviewed in the second section plays this role. However, due 
to the transfer dependence of the Sudakov factor, we must firstly face the problem of its 
analytic continuation [17[ . 

Before turning to this analytic continuation, we rewrite the expression of s from Eq. (|7|) 
for spacelike transfer, in the form: 

s(x,Q,b) = ^ln^L(U-l-lnU), U = (20) 
to explicitly see that s increases rapidly with U < 1 at large Q (remember that if x is 
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small, we can always consider the x-term in turn) so that the region of not too large 



suppression is U ~ 1, where we have: 



'"it^v-V- < 21) 



For large timelike transfer, setting Q = —iW in Eq. (EO), we get: 



1 4/3 16/? In 

with ssl and ?7 the spacelike expression of Eq. (EDI). In the above equation, the real 



s(x, W, b) - s SL (x, W, b)*i^-\nU- — ^ ^ (22) 



part is effectively small, while the imaginary part remains close to in the region of 
intermediate suppression. We can therefore assume, in the asymptotic regime, the same 
scale dependence for the Sudakov factor and use the expression ssl for our study. 

To simplify our purpose, we will concentrate on the simpler case one gets by consid- 
ering only the transverse behavior of the gluon propagator, that is setting the transverse 
momentum to k = in Eq. @. As we will show, this does not alter the naive behavior we 
previously get. However when taking the Fourier transform, only one transverse distance 
remains, b = hi = b 2 and after angular integration we are led to replace the integral I in 
Eq. ( |T7| ) by the quantity /': 

I' = [ A bdbK (-zb)e- 4s ^ b) (23) 
Jo 2 

where we have limited our study to the x = y = \ case. Thanks to the Sudakov suppres- 
sion, the upper bound of the integral is naturally b = A -1 . 

For the Sudakov exponent, we consider the approximate expression of section with 
the prescription of Li and Sterman [[/J which is to set the exponential to unity in the 
region that should not be controled by Sudakov evolution, here for b < l/\z\. With these 
simplifications, we get: 

/' = A (-1 - iK^-i) + [ W/2 uduK (-tu)e-^ w 'wA (24) 



i 

which is to be compared with the expression without Sudakov correction: 

J= P ( " 1_i T^ ( " w/2)) (25) 

We present in Fig. || the result of a numerical computation for both the real part (a) 
and the modulus (b) of the quantity 

W 2 

A = —I + l (26) 

which dictates the deviation from the counting rule canonical result and compare it to 
the original quantity ^-1 + 1. We observe that after an intermediate regime (W < 20A), 
the expressions including Sudakov suppression slowly decrease to contrary to the case 
of the rectangular wave function. 
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Another feature we have omited until now is the fact that the expression for the form 
factor is a superposition of amplitude with various fractions x, y weighted by smooth dis- 
tributions. This should also modify the W /_3//2 -power law and we examine this possibility 
in Appendix A. 

We conclude that the W~ 3 ^ 2 behavior does not resist to the inclusion of any realistic b 
or x-y integration procedure. Indeed the dimensional counting rules are valid at timelike 
transfers and furthermore the form factors are asymptotically the same. 



3.4 Comparison at intermediate energy 

Let us now turn to the discussion of the ratio of timelike over spacelike form factors 
in the intermediate range. The complete integration formula is: 

F = J dxdy J b 1 db 1 b 2 db 2 ^\x, 6i)^ (0) (y, h) T H e~ 5 (27) 



with the hard scattering amplitude from Eq. ([13]): Th = lQnasT. The integration range 
for the longitudinal fraction of momentum goes from to 1, whereas transverse distances 
it go from to 1/A thanks to the Sudakov suppression at large b. 

For the numerical study, we take into account the one loop running of the QCD 
coupling in the hard scattering: 

7T 

= (3 ln(t 2 /A 2 ) ' * = max (Vxy\z\,bx ,&2 ) (28) 

with the prescription for the renormalization point described in [|18| and A = 200MeV. 
The Sudakov factor e~ s contains the corrections for the two wave functions together with 
the anomalous running of the four quark operator Th [[18]. 

The Sudakov factor should be analytically continued as discussed in sub-section |3.3| [17] . 



It turns out that this procedure leads to quite mo del- dependent results at non asymptotic 
transfers. This has to do with the truncation of the exponentiated expression and with 
the need to suspect the validity of the approach when the (real part of the) Sudakov 
exponent becomes positive, i.e. when Sudakov suppression turns to an enhancement. 
We leave to Appendix B a somewhat detailed discussion of these effects, the conclusion 
being that an extra modification of the timelike value may come from this continuated 
Sudakov exponential, but that it is quite difficult to reliably quantify this statement. In 
the following, we will thus keep the Sudakov factor at its spacelike (real) value. 

We used various forms for the soft wave function ip^(x, b) to test the sensitivity of the 
result to this input. One may consider wave functions without intrisic transverse behavior 
(■^^(x, 6) = f(x)) and use either the asymptotic form: 




Vas{x) = ^-Ux{l-X), (29) 

with f w = 133MeV the pion decay constant, the CZ form: 
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<p C z(x) = 5(1 - 2xY<p aa (x), (30) 

or other expansions in terms of Gegenbauer polynomials like those proposed in ||, |5|. In 
the following, we will show some results for one form ||: 

ip FHZ (x) = (.6016 - 4.659(1 - 2xf + 15.52(1 - 2x) 4 )ip as (x) (31) 

As mentioned in section |2], the two last distributions have slow logarithmic evolution 
with Q. We ignore this evolution because it is quite insignificant in the range of energy 
we consider here. 

Following || [191, one may also include some intrinsic transverse behavior. We tried 



the different forms of wave functions described in [ylj and found only small differences 
for the behavior of the ratio. Therefore, we only quote here the sample form for which 
we will show some results in the following: 

^ (0) (x,6) = ip(x)e- b2/ < (32) 

which is a simple way to modelize the transverse behavior without a long set of parameters. 
2b related to the valence state radius is proposed in b% = 4.082 GeV~ 2 . 

Figs |7|-[8] show our numerical results for the meson form factors in the large but non 
asymptotic Q 2 = \q 2 \ region []. In Fig. [7], Q 2 \F n (Q 2 )\ is plotted against Q 2 for both timelike 
and spacelike regions up to Q = 50A. The distribution considered is the asymptotic one, 
Eq. The slow convergence of the timelike and spacelike quantities is manifest while 
the counting rule (F n oc l/Q 2 ) is reasonably well describing the Q 2 dependence down to 
a few GeV 2 in both cases. The inclusion of the intrinsic 6-dependence given by Eq. (|32"D 
(dashed lines) does not significantly modify the results. 

In Fig. [|(a), the modulus of the timelike form factor (multiplied by Q 2 ) is shown 
for the three choices of distribution amplitudes: CZ form (solid line), asymptotic one 
(dashed line) and FHZ one (long-dashed line). The experimental data shown comes from 
\I/ decay jnj. Sudakov suppression has been included but no intrinsic 6-dependence. 
Fig. |8|(b) shows the ratio of the timelike to the spacelike form factors. This ratio is rather 
wave function independent and decreases very slowly to 1 from a value of around 2 in 
most of the experimentally accessible range. 

Although this ratio turns out to be quite difficult to reliably extract from experimental 
data in the meson case, it is quite straightforwardly measured in the proton case. We will 
analyse the proton case in a forthcoming work. If we restrict to a simple quark-diquark 
picture, we would get a timelike to spacelike ratio quite similar to the one obtained here 
for the meson case, and thus understand the experimental value. 

^dn this sub-section, as there is no confusion, we do not distinguish the absolute value of spacelike and 
timelike transfers. 
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4 Conclusion 



In this paper, we showed that the difference between spacelike and timelike form fac- 
tors at large accessible transfer is predictible from an improved perturbative QCD analysis. 
We understand at least qualitatively the enhancement of the timelike values at large but 
non asymptotic transfers as mostly due to the integrable singularities of gluon and quark 
propagators. This strengthens the faith in the applicability of perturbative reasoning at 
intermediate energies (above a few GeV) at least for semi-quantitative understanding of 
the strong interaction physics. 

For the pion case, the uncertainties in the extraction of the spacelike form factor fT2|| 
show the need for another way to access this observable, the simplest one in exclusive 
scattering. We demonstrated that the usual formalism of Brodsky and Lepage has to 
be improved to account for the differences between timelike and spacelike regions in the 
energy range experimentally reachable. More experimental data are still needed to test 
our knowledge of the pion wave function. 

The proton case is more interesting since the extraction of the spacelike form factor is 
without ambiguity. The comparison of spacelike and timelike form factors thus appears 
to be a good way to understand the hadronic wave function. 

Many other hard exclusive processes dwell on timelike transfers. The 77 — > tctt, pp 
reactions at fixed angle for instance demand a more careful analysis than available now, 
not to speak of the difficult instances where pinch singular diagrams mix up, as in the 
ratio of pp to pp elastic scattering. More work needs to be done and experimentally 
tested before we know for sure that exclusive timelike reactions help us to understand 
confinement dynamics through the unraveling of hadron wave functions in their lowest 
Fock state. 
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Appendix A 



We give in this Appendix some arguments for the power suppression of the oscillating 
factor we get in sub-section |3.3| . Again we will concentrate on the simple case one has 
when neglecting the transverse momentum in the quark propagator (see sub-section [3.3|) 
and consider the integral analogous to I in Eq. (IT 



/" = [ B bdbKo(Jxyzb) = — ^ ^KAJxyzB). (33) 

jo xyz z wxyz 



As in section |3], we should worry about the behavior of the Bessel function with large 
argument ( Ki has the same asymptotic behavior as K ). We get: 



1 _ 

I" « - -^e~^ zB . (34) 

xyz 2 {Jxyzf 2 



In this explicit asymptotic form, one may guess that because the phase varies rapidly 
with x or y due to the presence of the large \z\ factor there may be some destructive inter- 
ferences when integrating the second term of I". The presence of any smooth distribution 
as weight functions will not destroy this feature. To see this explicitely we can perform 
the integration over some finite range for x to avoid region where the asymptotic form 
fails and also to allow further simplifications. Precisely, we look at: 

/•1/2+0 ff i 

J= dx<f(x)l"(x,y = -) (35) 

Jl/2-a 2 

and replaces x = 1/2 everywhere in the integrand except in the phase where we take the 
expansion of the square root of x around 1/2 up to first order. With these simplifications, 
we easily get: 



j 8a 1 N 
3 = T^2' 



7r _zBSinh^ 

— e 2 ^ 

zB a 



zBa " 

(36) 



which with the replacement z = —iW has a modified behavior compared to I" and the 
leading behavior is now identical in the spacelike and timelike directions of the complex 
plane. However, the previous equation, even if approximatively, still indicates qualita- 
tively that this identical asymptotic behavior may be reached rather slowly. 



Appendix B 

We discuss in this Appendix the analytic continuation of the Sudakov suppression 
factor. 
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In the complete expression for the form factor Eq. (p7|), two factors: T# and e~ s , are 
scale dependent and must a priori be analytically continued. In sub-section [T3| , we give 
some arguments to show that the Sudakov factor has a leading behavior which is not 
affected by analytic continuation so that we can study the contribution to timelike form 
factor ignoring this kind of difficulties. However, in the range of transfers we consider in 
the numerical study of sub-section |3j4], these arguments may not apply. In e~ s , the scale 



always appears in logarithms and after analytic continuation, one gets a phase which 
is sub-leading compared to the remaining large logarithm. The Sudakov exponent is 



known |15|] up to next-to-leading logarithms and we may keep the imaginary part which 
comes from the leading- log (Eq. (^)) as a next-to- leading component. In this kind of 
analysis, the additional real part due to the product of logarithms (the 7r 2 -factor), which 



may lead to an aditional enhancement ||17|| , is automatically dropped. A further analysis 
of the amount of correction which may be provided by such terms gives an 0(10%) extra 
enhancement of the timelike form factor. 

The effect of the imaginary part in the Sudakov exponent appears to be more impor- 
tant. The ratio of the timelike to spacelike factor is: 



e 



-SsL 



Cpir — In b 



^Wb)^-^^^. (37) 

We must define a prescription to cut-off the small b region. For the spacelike Sudakov 
suppression, Li and Sterman J7| choose to include the exponential factor only in the region 
of large b defined by bxW > \[2 and furthermore only when the sum of all exponents is 
negative and leads effectively to a suppression. The first prescription is associated with 
approximations discussed in section |2| and we will not questioned it in the following. In the 
spacelike case, it appears that the second constraint may be easily relaxed in a numerical 
study, because only a very small enhancement results when forgetting this constraint. We 
have observed that this is not likely to be the case for the timelike form factor. 

A numerical study, in the range of transfer 10-30A, with the consideration of the total 
gluon propagator alone and the prescription: from Eq. ([37]) if Re(— S) < 0, shows, for 
the asymptotic distribution, a 5-10% depletion of the timelike form factor with respect to 
the value without consideration of phase, whereas for the CZ distribution, the diminution 
is 20-30%. 
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Figure 1: The pion form factor. 



Figure 2: Factorization of the form factor. 



Figure 3: Tree graph for T#. 



Figure 4: Leading radiative corrections grouped in the wave function. 



Figure 5: The timelike hard amplitude. 



Figure 6: The real part (a) and modulus (b) of the deviation to spacelike scaling A. Solid 
(dashed) line is with (without) the Sudakov suppression factor. W is in A units. 



Figure 7: Timelike and spacelike form factors: different transverse behaviors. The x- 
dependence is the asymptotic form of Eq. (p9|). Energies are in A units. 



Figure 8: Timelike form factor (a) and timelike over spacelike ratio (b): sensitivity to 
distribution. CZ is Eq. (|30|), FHZ ([31]), as fl29|); experimental data from |TT . 
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